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We present a theoretical investigation of space separated electron 
and hole distributions, which consists in the self-consistent solving 
of the Schrodinger equations for electrons and holes and the Pois- 
son equation. The results are illustrated for the 

GaN/Alo.sGao.TN quantum well. The optical gain spectrum in 
a [0001]-oriented GaN/Alo.sGao.yN quantum well in the ultravio- 
let region is calculated. It is found that both the matrix elements 
of optical transitions from the heavy hole band and the optical 
gain spectrum have only the strict x (or y) light polarization. We 
present studies of the influence of the confinement of wave func- 
tions on the optical gain which implicitly depends on the built-in 
electric field calculated to be 2.3 MV/cm. Whereas the structures 
with narrow well widths exhibit the usual development of the light 
gain maximum almost without shifting the spectral region, a sig- 
nificant blueshift of the gain maximum is found with increase in 
the plasma density for wider quantum wells. This blueshift is as- 
cribed to the interplay between the screening of a strain-induced 
piezoelectric field and the bandstructure. A large Sommerfeld or 
Coulomb enhancement is present in the quantum well. 



1. Introduction 

Direct wide band gap group Ill-nitride semiconductors 
based on GaN and its alloys have received a great at- 
tention due to their applications in optoelectronic de- 
vices such as light-emitting diodes and lasers at green- 
blue and near-ultraviolet wavelengths and solar-blind 
photodetectors Pid- A number of ultraviolet light- 
emitting diodes (3|-|8| and laser diodes already 
have been demonstrated. Realizing the deep-ultraviolet 
semiconductor-based light-emitting diodes will provide 



compact high-efficiency light sources for various applica- 
tions, for example to the biological detection and data 
storage [^. Thus, these structures are in the develop- 
mental stage. 

Here, we present a theoretical investigation of the in- 
tricate interaction of the electron-hole plasma with a 
built-in electric field. For this purpose, the calculation 
of a quantum well bandstructure is performed using the 
invariant method and the envelope approximation. We 
consider a quantum well of width w in GaN, which is 
oriented perpendicularly to the growth direction (0001) 
and localized in the spatial region —w/2 < z < w/2. 
In the GaN/AlGaN quantum well structure, there is 
a strain-induced electric field. This piezoelectric field, 
which is perpendicular to the quantum well plane (i.e., 
in z direction) may be appreciable because of the large 
piezoelectric constants e in wiirtzite structures. 

The confinement of wave functions has a strong in- 
fluence on the optical gain which is observed with an 
implicit dependence on the built-in electric field which 
is calculated to be 2.3 MV/cm. Such fields are present 
in GaN/Alo.aGao.yN systems, because the strain is in- 
duced by the lattice mismatch. The relative magni- 
tude of piezoelectric effects depends sensitively on the 
quantum- well width and the plasma density. In this pa- 
per, we present the results of theoretical studies of the 
space separation of electron and hole distributions on the 
basis of the self-consistent solution of the Schrodinger 
equations for electrons and holes and the Poisson equa- 
tion. The Poisson equation contains the Hartree poten- 
tial which involves the space distributions of the charge 
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density for electrons and holes. We focus on details of 
the bandstructure for the sake of comparison of differ- 
ent quantum well structures. In the calculations of a 
band structure in the high-concentration regime, we dis- 
cuss the treatment of the quantum confined Stark ef- 
fect (QCSE). By comparing the gain spectra for two 
GaN/AlGaN quantum well structures with different well 
widths, we show the interaction of the bandstructure and 
the piezoelectric field. In particular, we will show that 
the wide-quantum-well structures, where the QCSE is 
appreciable, demonstrate a significant blueshift of the 
gain maximum, whereas the structures with a narrow 
well width exhibit the ordinary behavior of the light gain 
maximum almost without shifting the spectral region. 

A similar blueshift of the exciton resonance was 
observed and analyzed on the microscopic level for 
GalnN/GaN quantum- well systems They reflect 

a perturbation of the compensation between the self- 
energy and the field renormalization contribution to the 
microscopic interband polarization caused by a real spa- 
tial distribution of charges. Such a feature is charac- 
teristic of this quantum well and is not inherent to the 
GaAs quantum well die to the lacking of a piezoelectric 
field. Accounting the Coulomb renormalization of ma- 
trix elements of the electric dipole moment in the two- 
band model of quantum well structure causes a variation 
in the oscillator strength with a variation of the carrier 
density and the quantum well configuration. 

15fl , the matrix elements of the dipole mo- 



In work 



ment for interband transitions and the optical gain of 
a deformed wiirtzite GaN quantum well were presented 
without consideration of the intrinsic built-in piezoelec- 
tric field in the quantum well structure. 

In work the laser gain was investigated for Al- 
GaN wiirtzite quantum well structures. The optical gain 
spectrum was computed by simultaneously diagonalizing 
the kp Hamiltonian and by solving the Poisson equation. 
However, no significant shift of the gain maximum with 
increase in the plasma density in the framework of single 
structure was obtained. This indicates that, in the given 
structure in the high-density regime, QCSE is insignifi- 
cant. This result coincides with our calculations of the 
gain shown in Fig. 5. 

In work [ITj, a self-consistent calculation of the op- 
tical gain in pseudomorphically strained GaN quantum 
wells as a function of the carrier density was presented. 
But the spectrum renormalization and the electric dipole 
momentum which are caused by electron-electron and 
electron-hole Coulomb correlations were not considered 
there. 



Understanding the influence of the bandstructure and 
QCSE on laser gain properties should help one to im- 
prove the laser performance and the optimal configura- 
tions of a device. 

The light gain spectrum presented in the paper re- 
flects only the strict TE {x or y) light polarization. It 
is known that the valence-band spectrum at the 

r point originates from the sixfold degenerate Fis state. 
Under the action of the hexagonal crystal held and the 
spin-orbit interaction in wiirtzite crystals, Fis splits and 
leads to the formation of three spin degenerate levels: 
Fg, upper F7, and lower F7 levels. 

In Section 2 for the processes of emission or absorp- 
tion, we will calculate the energies and the wave func- 
tions of the lowest conduction subband and the valence 
subbands. The dependences of the matrix elements for 
dipole optical interband transitions and the light gain 
spectrum in GaN quantum wells on the quantum well 
width and the charge density are derived. Section 3 
presents the Hartree-Fock light gain spectra and the ma- 
trix elements for dipole optical interband transitions that 
are calculated within the theory described in Section 2. 
By comparing the light gain spectra for two GaN/AlGaN 
quantum well structures of different well widths, we show 
the interaction of the bandstructure, polarization field, 
and charge density. We determined the red renormaliza- 
tion of the light gain spectrum caused by the electron- 
electron and hole-hole Coulomb interactions. It is found 
that the Sommerfeld enhancement composes 26.7 gain 
value, which was obtained in the Hartree problem. This 
enhancement of the electric dipole momentum is caused 
by the electron-hole Coulomb attraction. 



2. Theory 

We consider QCSE in strained wiirtzite 
GaN/Alo.aCao.rN quantum wells with widths 2.6 
nm and 3.9 nm, in which the barrier height is a 
constant value for electrons and is equal to Uq = 490 
meV. The theoretical analysis of the optical gain of 
strained wiirtzite quantum well lasers is based on the 
self-consistent solution of the Schrodinger equations for 
electrons and holes in quantum well of width w with 
including Stark effect and the Poisson equation. The 
Poisson equation contains the Hartree potential which 
involves the charge density for electrons and holes. All 
researches are performed at a temperature of 300 K. 
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The first energy level of an electron in the quantum 
well of width w is equal to [23| 



El 



(1) 



where m = 0.19mo is an electron effective mass, and ^ 
is determined from equation 



cos^ = ±7^, 



(2) 



where 7 = tan^ > 0, and $ = For fco, the 

following equality holds: 



arcsm ■ 



(3) 



The wave function of an electron on the first energy level 
with regard for QCSE is as follows 2J|: 



VA 



(4) 



Here, 



M/(z,/3) = 
where ipi{z,j3) 



V-iCz,/?) 

V'(z,/3) 



(5) 



C sin {koz + So)e 



From the boundary conditions (23|j23 | i^i^z, P)\z=-w/2 



^J|g|.=»/2, we find Ci = Csin(-i^+<5o)e^** 



Sq = + arctan ^ , where A is the area of a quantum 
well in the xy plane, p is the two-dimensional vector 
in the xy plane, and fcj = (k^jky) is an in-plane wave 
vector. The constant multiplier C is found from the 
normalization condition 



\^{z,P)fdz = 1. 



(6) 



Such a representation of the wave function gives the in- 
formation that the conduction band corresponds to the 
Fy representation, which arises due to the splitting of 
the Cg^ space group by the crystal field with Fi. In 
other words, the conduction band wave functions orig- 
inate from S atomic orbitals. This is important at the 



derivation of matrix elements of the electric dipole mo- 
ment by the Wigner-Eckart theorem. 

The strong mismatch of the lattices in GaN and 
Alo.3Gao.7N leads to internal strains in the GaN layer. In 
noncentrosymmetric structures, the internal strains can 
induce a macroscopic built-in polarization field. This 
phenomenon is known as the piezoelectric effect. This 
phenomenon can also be described as a strain inducing 
an electric field. It is known that this piezoelectric field, 
which is perpendicular to the quantum well plane, can 
be significant because of the large piezoelectric constants 
in wiirtzite structures which are connected with one an- 
other: 



K \ \ C33/ 



sp ) J 



(7) 



where e is the piezotensor, Pgp is the spontaneous po- 
larization, e is the strain tensor, C13 and C33 are the 
elastic constants, and k is the permittivity of the host 
material. We calculated the built-in piezoelectric field 
in the GaN/Alo.3Gao.7N quantum well structure from 
relation (7) and found E ~ 0.23 x 10^ V/cm. 

We take [H, IHl the following values for con- 
stants: Ci3 = 106 GPa, C33 = 398 GPa, egi = 
-0.44 X 10^ V/cm, 633 = 0.66 x 10^ V/cm, P^p = 
—0.26 X 10^ V/cm. The transverse components of the 
biaxial strain are proportional to the difference between 
the lattice constants of materials of the well and the 
barrier and depend on the Al content: x, Cxx = ^yy — 

aAl3,Gai_^N — IGaN , n 

, aAUGai_^N — iGaN + X[aAlN — dQaN )', 

OGaN — 0.31892 nm, oain = 0.3112 nm. The longitu- 
dinal component of a deformation is expressed through 
elastic constants and the transverse component of a de- 
formation: Czz = -2§^ea;a;. 

One can find the functional, which is built from (4) 
and (5), in the form 



J(/?) = 



where 



(mm 



(8) 



(9) 



^2 ^2 

V{z) = J7(z) + e$(z). 



Hc = Eg +A1+A2 + ^^^t + acz<^zz + acJ^i^xx + ^yy), 

(10) 
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^ = 0.197710, and a^^^j^ = -4080 meV [30]. The 



quantity U (z) can be represented in the form 



e (-00... -w/2), 



U ^ ( eEz, z e [-w/2...u;/2], 

C/o + ^,ZG (7«/2...C30) 



(11) 



To account the piezoelectric effects, we modify the 
Schrodinger equation for electrons and holes, by in- 
cluding an off-diagonal contribution to the electron-hole 
Hamiltonian. The Schrodinger equation for an infinitely 
deep quantum well with regard for the QCSE and the 
Hartree potential created by spatially separated elec- 
trons and holes can be presented in the form 



(12) 



where iJ — Hj^ + eEz + e$(z). We introduce the Bloch 
function written as a vector in the three-dimensional 
Bloch space: 



\aav kt 



bi'\z,kt) 



|l,o-t,) 
\2,<Jv) , 

|3, (7^) 



(13) 



where 



(14) 



i=l 



and j — 1,2,3. The Bloch vector of the a- type hole with 
spin (T„ = ± and momentum kt is specified by its three 
coordinates [vl^^\n,a], V^^\n,a], V^'^\n,a]] in the ba- 
sis [|l,CTt,), |2,CTt,), |3, CT^)]. The envelope z-dependent 
part of the quantum well eigenfunctions can be deter- 
mined from the boundary conditions Xn{z = —w/2) = 
Xn{z = w/'i) = for an infinitely deep quantum well as 



/ 2 z 1 

Xn{z) = \/ - sin(7rn(- + -)), 
\ w w 2 



(15) 



where tt, is a natural number. The hole wave function 
can be written as 



*r"(r) 



gi kt p 

"7^ 



laav kt 



(16) 



where 7/ — {kt,a} in the envelope- wave approximation, 
in which the wave function is considered as a product of 
the envelope part x{z)e^^"' and a periodic Bloch multi- 
plier. The Bloch vectors in the envelope wave approxi- 
mation are projections of the exact Bloch vector on the 



subspace of vectors with the symmetry inherent to the 
r point [23. We have 



ff4 



F Kt T m 

Kt G A =F iHt 
±iHt A ± iHt A 



(17) 



in the basis [|1,(t^), |2,ct„), |3, cr.u)] [28'|, where 
F = Ai A2 + A + 61, 

G = Ai - A2 + A + 6*, 

A = Afe + Ae , 



9fc + i 



At 



-{Aikl + A2k'i), 



2mo 

Ae — DiCzz + D2{exx + ^yy): 



% = T^{A-iki + A^kt), 

% = D^ezz + Di{exx + ^yy), 

Ht = -^{Aektkz), 



Kt - 



2mo 



2mQ 
A = V2A3 



^ + ^y> 



|1,±) = i=[|l, 1) I t) e"^e-^ ± |1, -I) I i) e^e^], 
|2,±) = -i=[±|l, 1) I i) e^e-^ + |1, -1) I t) e^e^], 
|3,±) = ^[±|1,0) I t) e"^e-^ + |1,0) | ^e^^e^], 
\l,l) = -^\X + iY), 
|1,0)-|Z), 
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\\-\) = —\X-iY). 

The valence subband structure E^''{kt) can be deter- 
mined by solving the system of equations 



J2{H^i {K = ^i §-^) + SvE-'' (fc*)) fct) =0, (18) 



where i = 1,2, 3. In the quasicubic approximation, the 
parameters of effective mass and deformation potential 
are connected by the relation (l9l . 21 1 : 



4^5 - a/2A6 = A3, 2A4^-A3^Ai-A2, 



In calculations, we take the effective-mass parameters 
for the valence band [2^ as Ai — —6.56, A2 = —0.91, 
A3 = 5.65, A4, = -2.83, A5 = -3.13, Ag = -4.86, 
the parameters for deformation potential [30.j as Di = 
700 meV,D2 = 2100 meV, D3 = 1400 meV,D4 = 
—700 meV, and the energy parameters at 300 K |l5l . l25l | 



as E„ = 3507 meV, Ai 



A, 



16 meV, A2 



Aso/3 — 4 meV. Solving the Poisson equation 



47r 



(20) 



- V2D6 



2D4 



-D3 



D1-D2 



A2 A3. 



(19) 



with the condition p{z)dz = and with the selected 
wave functions, we find the Hartree potential e^{z): 



e<i> = 



2e2 



E 



^ ( COS TT Z{h-\-ri) 

w{- ^ 

M'2 



cos TT Z{rL—k) ^ 



2e' 



51 / ktdkt{S\{a,\C^\c7,)\S) hn{kt) 



-cos {-kQW+2So) ni e^'''0-'^)''+f ) 

2 4(ko-/3)^ 
-''g' 2cos2(fcoz+(5o)e~^'^"^ 



8/32 



(4/32+4feg)2 



1-cos {kow+2So) m g-2(^o+g)(»-f ) 
2 ^ 4(^H^^P 



sin 2(fcoz-f<5o)e ^'^ " ,^ o 



(21) 



r 



where Z — ^ + ^, ga and correspond to the de- 
generation of the a hole band and the first quantized 
conduction band, respectively, e is the value of electron 
charge, k is the permittivity of the host material, and 
fa,p{kt), /in(fct) are the Fermi-Dirac distributions for 
holes and electrons. Here we assume the charge concen- 
trations 9 X 10^^ cm^^, and 7 x 10^^ cm^^. 

Solving (12) for holes in the infinitely deep quantum 
well and finding the minimum of functional (8) for elec- 
trons in a quantum well with barriers of finite height, 
we can find the energy and the wave functions of elec- 
trons and holes with regard for the space distribution of 
electron and hole charge densities in the quantum well 
with given concentrations in a piezoelectric field. The 
screening field is determined by iterating Eqs. (8), (12), 



and (21) until the convergence of bandstructure calcula- 
tions is reached. We use the space carrier distribution of 
carriers in the lowest order for the envelopes of the wave 
functions of electrons and holes. 

Consider the matrix elements of interband transitions: 



(k) = J d^rU■i^„^kepUJcrk■ 



{22) 



The wave functions of the valence band transform ac- 
cording to the the representation Fi -f Fs, while the wave 
function of the conduction band transforms according to 
the representation Fi. In order to find the representa- 
tion for Mja-^ j'a-' {^), let us consider the direct product 
Fi X (Fi + F5). The symmetry elements of the point 
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group Cgt, are as follows: 
g = E, C2, 'i.C^i, 2C6, 3(Ti,, 3cr^, 



(23) 



where C„ is the axis of the n-th order, Scti, and 3cr(, are 
6 planes of reflection which pass through the sixth-order 
axis. For these elements, we find the representation Fi + 
F5: 



(24) 



(26) 



x(i?) = 3,x(C2) = -l,x(2C3)=0, 
x(2C6) = 2,x(3a„) = l,x(3cT;) = l. 

The squares of irreducible representation elements are 

g'' =E,E,C^,C:,,E,E. (25) 

We need to find 

X^(i?)=9,x|(C2) = l,x?,(2C3)=0, 
X^(2C6) =4,x|(3a„) = l,x|(3a;) = 1, 

whereas 

X^{E^) = 3,xv,(C^I) = 3,x^(2C|) = 0, 
X^(2C|) = 0,x^(3a2) ^ 3,xv.(3a;2) = 3. 

The symmetric representation can be found in the form 

\{xl{g) + x4g^))-- 
\{xl{E) + x^{E^)) = &, 
\{xl{C2) + xACl))^2, 

i(x^(2C3) + x^(2C|)) = 0, 
i(x|(2C6) + x42C62)) = 2, 
i(x|(3a,)+X^(3a2)) =2, 
Kx|(3<7;)+x43a;2)) = 2. 

The antisymmetric representations are 

\{xl{g)-x4g^))-- 
UxliE)-xAE'))-3, 



(28) 



Kx|(C2)-x^(c|)) 



^1, 



Mx|(2C3)-x^(2C|)) = 0, 



(29) 



i(x|(2C6)-xv.(2C|)) = 2, 
5(x|(3a.)-x^(3a2)) = -1, 
Kx^(3<)-x^(30) = -l. 

The symmetric representation can be decomposed into 
the irreducible representations 2 Ai + Ei + E2, whereas 
the antisymmetric one into A2+ Ei. Thus, the wiirtzite 



Hamiltonian -ff (e, k) must include the even functions 
(with respect to the time inversion), which are trans- 
formed according to 2Ai + Ei + E2, and odd functions, 
which are transformed according to A2 + Ei [l^ . 
The vector representation can be written as 



Xv{E) ^ 3,Xv{C2) - -1,X.(2C3) = 0, 
X.(2C6) = 2,x.(3a,) = l,x.(30 = 1, 



(30) 



and can be decomposed into the irreducible representa- 
tions Ai + El. The representation, according to which 
the interband operator is transformed, can be decom- 
posed into 



Fix(Fi+F5)=Ai+Si. 



(31) 



Thus, the direct product of representations 
(31) reflects the existence of nonzero matrix 
elements of the electric dipole moment of in- 
terband transitions, because the vector repre- 
sentation can be formed from these representa- 
tions. 

Allowed matrix elements of the electric dipole moment 
( 5|( crc|ep|wi, CTt,) are found in the form 



;5|(t |ep|fi,±) = -iP^e-^e-^Tsin^, 
;5|(t |ep|f2,±) = iP^e-^'^e^f sinfl, 
;5|(t|epb3,±) =±^F.e-*f COS0, 
[S\{i |ep|wi,±) = ±iPj_e^^'^e»T sin 61, 



(32) 



|ep|t-2,±) = T^P^e^^e-'^ sinO, 
(5|(;|ep|t-3,±) = -i=P,e'f COS0. 

Due to the symmetry properties of the Bloch functions, 
the only nonzero matrix elements between the basis func- 
tions are (2ll.[2^ 



;^|p,|i,o) = p, 

:^ip+ii,-i> = - 



[S\p-\1,1) = V2P^, 



(33) 



where ■p± = ± ipy Two constants of the matrix el- 
ements of the moment are as follows; P± = {S\px\X) 
and Pz = {S\pz\Z). Due to the cylindrical symme- 
try, the matrix element depends only on the difference 
<P = <y5E — <^k between the plane-projected angles of the 
vectors e||E and k. To simplify calculations, we assume 
(fik = and denote the spherical angles of the vector e 
by if and 9 [2l'|. We consider the case of a hole wave 
vector parallel to the c axis. In this situation, = in 
our calculations, and the vector e in the spherical coor- 
dinates takes the form e = (sin0cos(^,sin6'sin(/5,cos6'). 
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Fig. 1. Calculated square of the wave functions of a heavy hole (1) 
and a light hole (2) at the transverse wave vector fct = 8 X 10®cm~^ 
and an electron (3) (A); the charge density distribution on the 
quantum well 3.9 nm in width (B) at the charge concentration 



Fig. 2. Effective screening electric field distribution (C); the quan- 
tum well potential (1) and screening potential (2) on the quantum 
well 3.9 nm in width (D) at the charge concentration 9x 10^^ cm~^ 



whereas ep = isin0(e"^p_ + e~'^'^p+) +pz cos 9. It is 
known [21] that the values of constants |Pz,j_P can be 
found from the kp theory: 



, ^ 2 \{c\p. 
mo 



(34) 



From the value of experimentally measured conduction- 



band effective mass m 



2|p.._Lr 



15 eV. 



0.19mo and Eg = 3.5 eV, we 



obta 

■iit[) 

In Fig. 3, we show the fc-dependence of the matrix 
elements for the quantum well. We see that the matrix 
elements have the strict x (or y) light polarization for the 
transitions from the heavy hole band to the conduction 
band, while for the z light polarization, these transitions 
are forbidden [32]. That is why the light gain spectra 



presented in Figs. 4 and 5 reflect only the gain of TE 
polarized light for two widths of the quantum well. Such 
a behavior agrees with the results of calculations of the 
moment matrix elements for a wiirtzite GaN quantum 
well, since the valence band top originates from Fg, Fy, 
and r7 irreducible representations. The results which are 
presented in Figs. 3-5, testify to the optical polarization 
anisotropy of the matrix elements of the electric dipole 
moment for interband transitions in GaN/Alo.sGao.yN 
quantum well structures. 

The optical gain of a material 3 31 1 can be calculated 
from the Fermi golden rule 



ao 



IT e 



'Kmo wu) 
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Fig. 3. Moment matrix elements for the x- (or y-) polarization (A) 
and the z-polarization (B): (la) quantum well 3.9 nm in width at 



(2a) quantum well of 
p = 7 X 10^2 cm-2; (16) 



width 3.9 nm at the concentration n 
quantum well 2.6 nm in width at the concentration n = p = 9 X 
10^^ cm~^; (26) quantum well 2.6 nm in width at the concentration 
n = p = 7 X 10^2 cm-2 
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where e is the magnitude of electron charge, mo is the 
electron rest mass in the free space, c is the velocity of 
light in the free space, k = 8.27 is the permittivity of 
the host material, and are the Fermi-Dirac dis- 
tributions for electrons in the conduction and valence 
bands, respectively, e is a unit vector of the vector po- 
tential of the electromagnetic field, ma(^t) the in- 
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Fig. 4. Calculated Hartree gain spectrum: (la) quantum well 
3.9 nm in width at the concentration n = p = 9 X 10^^ cm~-^; (2o) 
quantum well 3.9 nm in width at the concentration n = p = 7 x 

lO^^cm^^; (16) quantum well 2.6 nm in width at the concentration 
(26) quantum well 2.6 nm in width at the 
■-p = 7 X 10^2 cm-2 
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Fig. 5. Calculated Hartree-Fock gain spectrum for the quantum 
well 2.6 nm in width and at concentrations: n = p = 7x 10^^ cm~^ 
(i) and n = p = 9 X lO^^ (,^-2 ^gj 



terband energy of the conduction and valence bands, hoj 
is the optical energy, and 7 is a half-linewidth of the 
Lorcntzian function, which is equal 6.56 meV. We con- 
sider the electromagnetic wave which propagates in the 
plane of the quantum well. The modal gain, which deter- 
mines the threshold condition of a laser, is proportional 
to the material gain multiplied by the optical confine- 
ment factor r and by the number a of quantum wells in 
the case of multiple quantum wells: a — aoT a. We take 
r equal to 0.01, and a is taken to be 1 in the calcula- 
tions. 
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Although the carriers within each band are in a 
strongly nonequilibrium state, the interband relaxation 
times are much larger than intraband relaxation times. 
Therefore, the Fermi-Dirac statistics can be used in the 
calculations. 

Using the expressions for the basis functions, we ob- 
tain two scalar polarizations for the matrix elements of 
the electric dipole moment. For the TE-poIarization 
(e = a; or y_L c axis), i.e., for the light polarization vector 
lying in the quantum well plane, we have 

|(M.)?Jfc,)P = 

= l<^l^-l^)l^ {(^i(/3)|E„V-,i[n,a]xn)^ + 

for£7 = +, (36) 

+(*i(/3)IEn^/f,[ri,a]x„)n, 
fora = — . 

For the TM-polarization (e = 5||c axis), i.e., for the light 
polarization vector, which is perpendicular to the quan- 
tum well plane, we have 

forcr = +, (37) 
for cr = — . 

3. Results and Their Discussion 

To describe the interplay of the bandstructure and the 
polarization effects in the Hartree problem, we consider 
the 2.6-nm and 3.9-nm GaN/Alo.sGao.rN quantum well 
structures. In the quantum well 2.6 nm in width at a 
concentration of 9 x 10^^ cm~^, a optical gain maximum 
is equal to 23673.7 cm^^ at the wavelength A — 334.5 
nm; while, at a concentration of 7 x 10^^ cm~^, the op- 
tical gain maximum is equal to 14245.9 cm^^. Such a 
gain is observed at the wavelength A = 336 nm. In 
the quantum well 3.9 nm in width at a concentration of 
9 X 10^^ cm^^, the optical gain maximum is equal to 
5752. 7cm^^ at the wavelength A = 348.3 nm; while, at 
a concentration of 7 x 10^^ cm~^, the optical gain maxi- 
mum is equal to 2840.8 cm^^. Such a gain is calculated 
at the wavelength A = 352.1 nm. Thus, the optical gain 
in the GaN/Alo.aGao.rN quantum well develops in the 
ultraviolet spectral region, as shown in Fig. 4. 



Numerically solving the Schrodinger equations (8) and 
(12) for electrons and holes and the Poisson equation 
(20), the steady state solutions allow us to construct 
the squares of the wave functions of a heavy hole ( i ) , 
a light hole {2) (e.g., at the transverse wave vector 
fct = 8 X 10^ cm~^), and an electron (5) (A); the charge 
density distribution on the quantum well width (B); the 
effective screening electric field distribution (C); and the 
quantum well potential and the screening potential on 
the quantum well width (D). The results of calculations 
for the concentration 9 x 10^^ cm^^ for a 3.9-nm quan- 
tum well are shown in Figs. 1 and 2, and the Hartree 
gain spectra are presented in Fig. 4. For the narrow 
quantum well, we see that, at a density of 7 x lO^^cm"^, 
the light gain is gradually developed, as the carrier den- 
sity increases. At high densities (i.e., when the density is 
equal to 9 x 10^^ cm~^), the optical gain develops nearly 
in the spectral region of the original optical gain at a 
plasma density of 7 x 10^^ cm~^. 

The behavior of the light gain coefficient for two quan- 
tum well widths and at given concentrations can be un- 
derstood from Figs. 1(A,B) and 2(C,D). From Fig. 1(A), 
we see that the overlap between the quantum confined 
electron and hole wave functions is related to the charge 
density distribution over the quantum well width, which 
is shown in Fig. 1(B). We can conclude that, for the wide 
quantum well, the overlapping integral is smaller than 
that for a narrow quantum well and reduces stronger 
with decrease in the carrier density. 

The effective screening electric field distribution 
for wide-bandgap GaN/AlGaN quantum well systems, 
which is presented in Fig. 2(C), is similar to that of the 
electric field in a condenser. 

As shown in Fig. 4, the situation is quite different for 
the 3.9-nm GaN/Alo.aGao.rN quantum well structure. 
Because of a weaker quantum confinement in this rela- 
tively wide quantum well, the piezoelectric field is able 
to significantly reduce the overlap between the quantum 
confined electron and hole wave functions, which can 
be seen from the comparison of Figs. 1 and 2. As a 
result the interband dipole matrix clement or the oscil- 
lator strength is substantially smaller than that in the 
case for the narrow 2.6-nni quantum well. This intrinsic 
quantum confined Stark effect can also cause a signifi- 
cant rcdshift of the gain maximum at a plasma density 
of 7 X 10^^ cm~^ as compare with the flat-bottom band 
situation. As the plasma density increases, the screening 
of the QCSE increases the overlap of electron-hole wave 
functions and, hence, the exciton oscillator strength. Si- 
multaneously, the weakened piezoelectric field, which in- 
duced earlier the redshift, leads to the net of blueshifts 
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in the gain maximum and the absorption edge with in- 
crease in the plasma density, as shown in Fig. 4. 

To calculate the concentration dependence of many- 
body Coulomb effects in the absorption spectrum of a 
GaN quantum well, we apply the method developed in 
[3, [33, [34 1 . Numerically solving the microscopic polar- 
ization equation, we see that, with increase in the plasma 
density, the optical gain (i.e., negative absorption) de- 
velops in the spectral region of the original exciton reso- 
nance. With increase in the free-carrier density, the ion- 
ization continuum shifts rapidly to the long-wavelength 
side, whereas the Is-exciton absorption line stays al- 
most constant, due to a high degree of compensation 
between the weakening of the electron-hole binding en- 
ergy and the band-gap reduction, like that earlier found 
for GaAs [35]. At high electron- hole concentrations, the 
electric dipole moment renormalization effects give rise 
to a large optical gain which is shown in Figs. 5 and 
6. The maximum of the Hartree-Fock gain spectrum 
equal to 29991.2 cm~^ is observed at the wavelength 
A — 339.1 nm at a concentration of 9 x 10^^ cm~^. 
From Fig. 6, we see that the Hartree-Fock spectrum 
is shifted to the long-wavelength side relative to the 
Hartree gain spectrum. Moreover, a large Sommelfeld or 
Coulomb enhancement is present in the quantum well. 
It is caused by an increase of the oscillator strength due 
to the electron-hole Coulomb attraction. 

The exchange Hartree-Fock energy spectrum renor- 
malization is accounted in the equation of motion for 
the microscopic dipole of the electron-hole pair. For high 
concentrations, this value is significant. It is somewhat 
larger for electrons and less for holes. Totally, this value 
is reflected in the Hartree-Fock gain shifting in compar- 
ison with the Hartree spectrum, as shown in Fig. 6. It 
should be noted that the gain spectrum involves not only 
Hartree-Fock correlations, but correlations of higher or- 
ders in the expansion in the Coulomb potential energy. 
This is achieved by the summation of the series in the 
Coulomb energy in the microscopic polarization equation 
in all orders of perturbation theory. In more details, the 
microscopic polarization equation for the dipole of an 
electron-hole pair for the wiirtzite quantum well will be 
considered in our next paper. 

4. Conclusions 

In summary, the self-consistant calculations of the 
Schrodinger equations and the Poisson equation of wide 
bandgap GaN/AlGaN quantum well systems show the 
interesting dependences of the matrix elements for dipole 
optical interband transitions and the light gain spectrum 
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Fig. 6. Calculated Hartree-Fock gain spectrum for the quantum 
well 2.6 nm in width at the concentration n = p = 9 X lO'^^ cm~^ 
(1) and the Hartree spectrum at the concentration n = p = 9 X 
10l2 cm-2 (2) 



on the quantum well width and the charge density. A 
blueshift with increase in the plasma density in the gain 
spectrum in relatively wide wells occurs as a consequence 
of the screening of the piezoelectric field induced by the 
quantum confined Stark effect, whereas the structures 
with narrow well widths exhibit the usual dependence of 
the development of the light gain maximum almost with- 
out shifting the spectral region. It is found that the ma- 
trix elements of optical transitions from the heavy hole 
band have the strict TE light polarization like that of 
the light gain spectrum. With regard for the Coulomb 
interactioin, a red shift of the Hartree-Fock light gain 
spectrum relative to the Hartree gain spectrum and a 
large Sommerfeld enhancement for a quantum well are 
found. 

The author is grateful to Prof. V.I. Sheka and Prof. 
V.A. Kochelap for numerous discussions. 
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OnXHHHA nOJIJIPHSAIlIHHA AHI30TPOni51, 

BHYTPimHIH E<1>EKT UITAPKA KBAHTOBOFO 
KOHOAHHMEHTy I BnjIHB KYJlOHIBCbKHX 
EOEKTIB HA J1A3EPHI XAPAKTEPHCTHKH 
[0001]-OPieHTOBAHHX GaN/Alo.aGao.rN 
KBAHTOBHX 5IM 

JI.O. JIoKomb 

P e 3 lo M e 

y ij,iH CTaTTi npe^CTaBjieno xeopeTH'^He flOcjiifl^KeHHH npocTopoBO 
po3flijieHHx ejieKTpoHHHx i flipKOBHx posnoflijiiB, aKe Biflo5pa>Ka- 
eTbca y caMoysro/i^eHOMy posB'HsaiiHi pisHaHb lIlpefllHrepa fljia 
ejieKTpoHlB Ta flipoK 1 piBHaHHa IlyaccoHa. PesyjibTara npoijiiocT- 
posaHO fljia GaN/Alo sGao^yN KBaHTOBoi smh. CneKTp onTHHHo- 
ro niflCHjieHHH b [OOOlJ-opieHTOBaHoi' GaN/AIo,3Gao,7N KBaHTO- 
Boi' aMH o6MiicjieHO b yjibTpa4)iojieTOBiii o6jiacTi. 3naii^eHO, mo 
aK MaTpHHHi ejieMBHTH onTHMHHx nepexofliB 3 Ba>KK01 /^ipKOBOl 
niflsoHH B 30Hy npoBiflHOCTi, TaK 1 cneKTp onTH^iHoro niflCHjieHHa 
MaioTb cTporo X (a6o y) nojiapH3ax(iK) cBiTjia. HoKasaHO ehjihe 
KOH4)aHHMeHTy XBHjibOBHX (J)yHKi;ift Ha onTHHHe niflCHjieHHa, axe 
HeaBHO 3ajie?KHTi> bI^ B6y^OBaHoro ejieKTpHMHoro nojia, mo o6mh- 
cjieHe i flopiBHioe 2,3 MB/cm. 51kii],o cxpyKTypH 3 By3i>KOio mnpn- 

HOio aMH npoaBJiaiOTb 3BHHaHHy 3ajlC>KHiCTb p03BHTKy MaKCHMy- 

My niflCHjieHHa CBlTjia MaiijKe 6e3 3MimeHHa cneKTpajibHoi' o5jia- 
ctI, to 3HaHHoro rojiy6oro 3MiiLi;eHHa MaKCHMymy niflCHjieHHa 3i 
3pocTaHHaM rycTHHH njiaaMH Ha6yBaioTb cTpyKTypn 3i 3HaHHoio 
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niHpHHOio KBaHTOBoi' 5IMH. LI,e rojiySe 3Miiii,eHHH Bi^HOCsTt B3ae- 4)opMaii,ieio i sohhoio CTpyKTypoio. BejiHKe 30MMep4)ejiBfliBCi>Ke 
MOflii MijK eKpaHyioHHM n'e3oejieKTpHHHHM nojieM, cTBopeHHM fle- a6o KyjioniBCbKe niflCHjieHHa npncyTHe y KBaHTOBiii aMi. 
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